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Abstract 


Previous versions of GISS climate models have either used formulations of 
Rayleigh drag to represent unresolved gravity wave interactions with the model resolved 
flow or have included a rather complicated treatment of unresolved gravity waves that, 
while being climate interactive, involved the specification of a relatively large number of 
parameters that were not well constrained by observations and also was computationally 
very expensive. Here, we introduce a relatively simple and computationally efficient 
specification of unresolved orographic and non-orographic gravity waves and their 
interaction with the resolved flow. We show comparisons of the GISS model winds and 
temperatures with no gravity wave parametrization; with only orographic gravity wave 
parameterization; and with both orographic and non-orographic gravity wave 
parameterizations to illustrate how the zonal mean winds and temperatures converge 
toward observations. We also show that the specifications of orographic and non- 
orographic gravity waves must be different in the Northern and Southern Hemispheres. 
We then show results where the non-orographic gravity wave sources are specified to 
represent sources from convection in the Intertropical Convergence Zone and 
spontaneous emission from jet imbalances. Finally, we suggest a strategy to include 
these effects in a climate dependent manner. 
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1 . 


Introduction 


The GISS Model-E climate model has been described in Schmidt et al. (2006), 
and this was the GISS model that generated results that were used in IPCC (2007), the 
most recent IPCC assessment. There were three versions of GISS Model E that were 
described in Schmidt et al. (2006). Two of these, M20 and F20 had model tops at 0.1 hPa 
and had 20 layers in the vertical, with M20 having 4 0 (latitude) x 5 0 (longitude) 
horizontal resolution and F20 having 2 0 x 2.5 0 horizontal resolution. The other, M23, 
had 4 0 x 5 0 horizontal resolution and 23 layers, with a top at 0.002 hPa. The Arakawa 
B-grid was used in all these models, and a sigma vertical-coordinate was used up to 1 50 
hPa, with a pressure vertical coordinate above. The basic model physics for these models 
were described in Schmidt et al. (2006) and previous GISS publications referenced 
therein. Schmidt et al. (2006) described how the results lfom these models compared to a 
variety of diagnostics from observations. 

Although most of the model physics used in these models were pretty much state- 
of-the-art, the treatments of unresolved gravity waves were not. Both the M20 and F20 
models used a Rayleigh drag scheme at the model top together with another simple 
Rayleigh drag scheme in the model interiors, above 1 50 hPa. While, these 
parameterizations were meant to crudely represent the influence of unresolved gravity 
waves, their treatments were not self-consistent in this regard. The M23 model used a 
climate-dependent gravity wave drag due to Rind et al. (1988). This included treatments 
of orographic drag and penetrating convection, shear and deformation gravity wave 
sources. While this latter gravity wave treatment adjusts to different climate regimes, it is 
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expensive in terms of computer time and contains a large number of adjustable 
coefficients that have little in the way of observational constraints. 

GISS is now readying themselves to participate in IPCC AR5. This has motivated 
a lot of updating of the various physics packages, and we have been implementing a new 
gravity wave treatment in a new GISS climate model. In the following, we will describe 
this, as well as show some comparisons with observations, and discuss how these new 
gravity wave treatments improve upon the Rayleigh drag parameterization previously 
used in GISS Model-E. 

2. Gravity Wave Effects 

Since the pioneering work of Leovy (1964) and the subsequent papers of 
Schoeberl and Strobel (1978) and Holton and Wehrbein (1980), the importance of 
parameterizing the effects of unresolved gravity waves to successfully model the middle 
atmosphere circulation has been apppreciated. Later, Palmer et al. (1986) and McFarlane 
(1987) showed that the effects of unresolved, orographically forced gravity waves should 
be parameterized in order to obtain good simulations of the troposphere. Even with the 
inclusion of orographic gravity wave parameterizations, however, it is still necessary to 
parameterize the effects of gravity waves arising ffom non-orographic sources, such as 
convection, frontogenesis, and jet sources in order to obtain reasonable middle 
atmosphere climatologies in climate models with reasonable resolution. There exist high- 
resolution atmospheric general circulation models that give realistic atmospheric 
structure without any treatment of unresolved gravity waves ( e. g., Watanabe et al., 
2008), but these models are too expensive to run with extensive climate interactions (e. 
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g., ocean, cryosphere, biosphere) using the present generation of computers. Thus, 
climate models will require parameterizations of unresolved gravity waves for at least the 
next decade or so. 

The development of gravity wave parameterizations began with the classic work 
of Lindzen (1981), and this was followed by several different formulations for 
parameterizing non-orographic gravity waves (e. g., Hines, 1997; Alexander and 
Dunkerton, 1999; and Warner and McIntyre, 2001). There have been more recent 
formulations for orographic gravity wave parameterizations (e. g., Lott and Miller, 1997, 
and Scinocca and McFarlane, 2000), and there have also been recent efforts toward 
including physically-based, non-orographic gravity wave treatments (i. e., that use the 
modeled phenomena in the climate model for sources of gravity waves in their 
parameterization) in papers such as those by Charron and Manzini (2002) and Richter et 
al. (2010). 

Our efforts have been motivated by the following considerations. We want 

1 . methods for gravity wave parameterization that give realistic 
atmospheric structures, interannual variability, responses to climate 
perturbations, and realistic transport characteristics. 

2. parameterizations that are physically well-founded. 

3. parameterizations that are computationally efficient. 

4. parameterizations that can be adapted to respond to a changing 
climate. 


5 



The previously used J-drag in GISS Model-E (Schmidt et al., 2006) does not 
satisfy the above criteria in that, although it is meant to simulate gravity wave effects, it 
behaves quite differently physically. The J-drag in Model-E responds to local winds 
rather than to wind filtering by levels below. Also, the J-drag provides a relaxation 
toward a motionless state for both the zonally averaged wind and for the zonally 
asymmetric wind. The climate responses to models with Rayleigh drag have been 
investigated by Shepherd et al. (1996), as well as by Shepherd and Shaw (2004), and they 
found that models with Rayleigh drag gave spurious climate responses. The desirability 
of momentum conservation for gravity wave parameterization was investigated in a series 
of papers by Shepherd and Shaw (2004), Shaw and Shepherd (2007), and Shaw et al. 
(2009). While it is true that Model-E did seek to conserve momentum by balancing the 
momentum deposition above with a counterbalancing momentum deposition below, the 
manner in which this is done is rather arbitrary, and does not physically correspond to the 
workings of gravity waves. 

Our initial efforts use the orographic gravity wave parameterization of McFarlane 
(1987) and the non-orographic scheme of Alexander and Dunkerton (1999). These 
choices were motivated by the following considerations. While the Lott and Miller 
(1997) and Scinocca and McFarlane (2000) schemes are more realistic than that of 
McFarlane (1987) in that they include such effects as low level wave breaking, upstream 
blocking, and lee- vortex dynamics, they are also more complicated to code for use in the 
GISS models, whereas the McFarlane (1987) scheme involves only simple coding. 
Furthermore, while some improvements have been noted when these newer schemes are 
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used, these improvements are relatively modest compared with the improvements that are 
realized when our gravity wave parameterizations are used (e. g., Scinocca et al., 2008) 
instead of the earlier Rayleigh drag treatments in GISS Model-E. 


Our use of the Alexander and Dunkerton (1999) gravity wave scheme was 
motivated by results in McLandress and Scinocca (2005). They showed that differences 
in non-orographic gravity wave parameterization schemes were less important than were 
the proper specifications of the nature of gravity wave sources (e. g., the gravity wave 
source spectrum). This, together with the fact that the Alexander and Dunkerton (1999) 
gravity wave parameterization allows for particularly simple mapping of the nature of the 
source spectrum to their effects on the mean flow, motivated our choice for the non- 
orographic parameterization. 


3. Some, simple, intermediate results 

There are several parameters that need to be specified for even our relatively 
simple choice of schemes. The McFarlane scheme uses variances of elevation calculated 
from a high-resolution topography data set, but one must specify values for two 
parameters. One of these is the critical Froude number Fr c , which determines the 
threshold for gravity wave breaking, and hence the vertical distribution of the gravity 
wave momentum deposition . We have taken Fr c to be 0.5, a conventional value. The 


• • • // 2 
other, in McFarlane ’s notation, is E-^-h e , 


where p e is a characteristic horizontal 


wavenumber and h e is a characteristic wave amplitude, and E is a constant that is meant 
to represent the wave intermittency. For the results shown in this paper, we have taken 
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E^- to be equal to 5.5 x 10' 6 m' 1 , and the h e 2 is taken from the topography height 


variances in each grid box. 


In the Alexander and Dunkerton scheme, we must specify the shape of the gravity 
wave spectrum, and this includes specification of a functional form as well as a width 
parameter. We are using the B 2 spectral shape of Gong et al. (2008) with the width 
parameter Cw equal to 10 m s' 1 and the source amplitude B m = 0.01 m 2 s' 2 everywhere in 
this paper. We use four azimuthal directions (North, East, South, and West) in which the 
gravity waves are launched. Finally, we launch our non-orographic waves at 100 hPa. 
This choice is motivated by two considerations. One is that jet imbalances are known to 
be a source of gravity waves (see Gong and Geller, 2010, for example) and the other is 
that deep convective towers impinging on the tropopause are also known to be a 
significant wave source, particularly in the tropics. 

In practice, there is a great deal of “tuning” that goes into the choice of gravity 

wave parameterization parameters since until recently there have been few observations 

to guide their choice, although this situation is now changing (see Alexander et al., 2010) 

with the growing literature on observations of gravity waves by various techniques. Both 

* For consistency and clarity, we use the terminology of Fr as in McFarlane (1987), which is actually the 
inverse Froude number (Scinocca and McFarlane, 2000). 
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resolved waves and unresolved gravity waves influence the atmospheric zonal mean wind 
and temperature states. Given that there are several adjustable parameters in both the 
McFarlane (1987) orographic gravity wave scheme and in the Alexander and Dunkerton 
(1999) gravity wave scheme, our philosophy in choosing values for these parameters is to 
first “tune” the orographic scheme to get the troposphere/lower stratosphere to agree 
reasonably with Northern Hemisphere winter zonal mean temperature and wind 
observations and then to “tune” the Alexander and Dunkerton (1999) non-orographic 
gravity wave scheme to agree with wind and temperature observations in the upper 
stratosphere. 

Our discussion in this section is based on a succession of four figures, each 
comparing ERA-40 with GISS model results that include no specification of gravity wave 
drag (referred to as GISS-ND); GISS model results with the effects of orographic gravity 
waves only included (referred to as GISS-OG); and finally GISS model results including 
both orographic and non-orographic gravity wave effects (GISS OG&NOG). These 
simulations have been carried out in an AMIP (Atmospheric Model Intercomparison 
Project) sense. We have run the model starting in 1979 for 21 years, and we compare our 
model climatology for the years 1980-1999 with that from ERA-40 for those same years. 
The horizontal resolution for all the GISS model results shown here is 2° x 2.5° as in 
GISS Model-E F20 in Schmidt et al. (2006), and the model tops are at 0.1 hPa with 40 
layers in the vertical. In these figures, the momentum flux, B t , at the 100 hPa source 
level was specified as 0.0015 kg m' 1 s' 2 at each azimuth. Sea-surface temperatures, ice 
conditions, ozone, and greenhouse gas concentrations are specified for the modeled years. 
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Figure (1) shows results for January zonal mean zonal winds, figure (2) results for 
January zonally-averaged temperatures, figure (3) for July zonal mean zonal winds; and 
figure (4) for July zonally-averaged temperatures. 

Looking at figure (1), note that the GISS-ND January subtropical jets in both 
hemispheres compare well with the ERA -40 results in both speeds and latitudes, and the 
easterly summer jet also looks somewhat reasonable although there is no evidence in 
ERA-40 of the double jet structure seen in GISS-ND above about 5 hPa. The GISS-ND 
winter westerly jet is much too strong, however, with zonal mean zonal winds of about 
100 m s' 1 above at 1 hPa, while ERA-40 only has winds of about 45 m s' 1 . Note also that 
the region of relatively weak westerlies between the winter upper and lower jet structures 
are always above 25 m s' 1 in GISS-ND while the ERA-40 winds are less than 20 m s' 1 in 
this region. The inclusion of orographic drag improves the agreement between our 
modeled January zonal mean winds and observations. The maximum westerly winds at 1 
hPa in GISS-OG are now about 50 m s' 1 , which is only a bit greater than the ERA-40 
winds in this region. Note also that the 30 m s' 1 , contour is at about 25 hPa and at about 
65 °N, which agrees well with ERA-40, whereas in GISS-ND, the 30 m s' 1 contour was at 
about 50 hPa and at about 60 °N, so the inclusion of orographic gravity wave drag has 
reduced the shear in the region between the jets to agree well with observations. With the 
non-orographic drag included, the maximum westerly winds at 1 hPa are about 45 m s' 1 , 
and they are located at about 40 °N. This is a bit equatorward of what is seen in ERA-40. 
In GISS OG&NOG, the minimum winds between the tropopause and polar night jets are 
below 20 m s' 1 , again in agreement with ERA-40. Looking at January summer, however, 
we see that the maximum easterly winds at 1 hPa are about 45 m s' 1 , which is less than 
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the ERA-40 observed 65 m s' 1 , but their latitude agrees well with observations. 

In figure (2), consistent with the thermal wind relation, the January winter lower 
stratosphere temperatures are too cold in GISS-ND, being below -90 °C 
compared to observed values of below -70 °C. This winter cold bias extends upward 
through the stratosphere, and there is also a summer warm bias in the GISS-ND so 
that the 1 hPa pole-to-pole temperature gradient is about 85 °C in GISS-ND 
and only about 40 °C in ERA-40. Consistent with the thermal wind relation, the inclusion 
of orographic drag has raised the minimum January winter polar night temperatures from 
-90 °C in GISS-ND to about -80 °C, which is closer to the ERA-40 polar night 
temperatures of -70 °C. It also has lowered the pressure altitude of the winter polar 
temperature minimum from about 20 hPa in GISS-ND case to about 40 hPa in GISS-OG, 
which is closer to the observations (~50 hPa). In the January winter hemisphere, the 
minimum polar night temperatures in GISS OG&NOG are about -70 °C, which are close 
to those in ERA-40, and their pressure altitudes compare well. At 1 hPa, the pole-to-pole 
temperature gradient is about 40 °C, which compares well with ERA-40. 

Looking at figure (3), the July GISS-ND subtropical jet latitudes and wind speeds 
compare well with ERA-40. At higher altitudes, both the summer and winter wind 
speeds are too high in GISS-ND, with the maximum modeled winter westerlies being 
about 175 m s' 1 at 1 hPa compared to the 95 m s' 1 maximum ERA-40 wind speeds at 1 Pa. 
The summer easterlies in GISS-ND are about 65 m s' 1 at 1 hPa whereas the ERA-40 
winds there are about 40 m s' 1 . The July winter westerlies in GISS-OG are decreased a 
bit. The maximum westerly winds at 1 hPa in GISS-ND are about 175 m s' 1 , and are at 
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about 65 °S. In GISS-OG, the maximum July westerlies at 1 hPa are about 150 m s' 1 , and 
are approximately at 60 °S. This is compared to westerlies of about 95 m s' 1 at 1 hPa in 
ERA-40 that are at about 45 °S. Also, note that the minimum winds between the 
tropopause jet and the polar night jet in both GISS-ND and GISS-OG are in excess of 30 
m s' 1 whereas in ERA-40, they are less than 30 m s' 1 . The GISS-OG&NOG July 
maximum winter westerlies are about 135 ms' 1 , which is about 40 m s' 1 more than in 
ERA-40. Also, these maximum westerly winds are at about 60 °S, compared to about 50 
°S in ERA-40. The observed equatorward tilt of the polar night jet seen in ERA-40 is not 
evident in the GISS-OG&NOG results shown in figure (3). The minimum winds 
between the tropopause jet and the polar night jet are above 30 m s' 1 in GISS-OG&NOG, 
which is more than in ERA-40. Looking at the July summer easterlies in GISS- 
OG&NOG, we see that the maximum winds are 55 m s' 1 and are located at about 15 °N. 
This is to be compared to the ERA-40 value of about 40 m s' 1 , which is located at about 
30 °N. 


Figure (4) shows the temperature comparisons for July. The coldest winter lower 
stratospheric temperatures in GISS-ND are about -1 15 °C whereas the corresponding 
ERA-40 temperatures are about -90 °C, and the pole-to-pole temperature gradient at 1 hPa 
in GISS-ND is more than 130 °C, whereas the corresponding ERA-40 temperature 
gradient is again about 40 °C. The inclusion of orographic gravity wave effects in GISS- 
OG has greatly improved the lower stratosphere polar night temperatures, with minimum 
stratospheric temperatures of about -95 °C centered at about 20 hPa in GISS-OG. The 
ERA-40 minimum temperatures are about -90 °C and are at about 30 hPa. The minimum 
July stratospheric polar night temperatures in GISS-OG&NOG are about -95 °C and are 
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at about 20 hPa. This is a little higher and colder than in ERA-40. The pole-to-pole 
temperature gradient at 1 hPa in GISS-OG&NOG is about 70 °C, which is about 30 °C 
more than in ERA-40, and is consistent with the stronger winds than in ERA-40. 

Summarizing the results shown so far then, it is apparent that the inclusion of 
orographic gravity wave effects has brought both the January and July zonal mean zonal 
winds and temperatures into much closer agreement with observations, particularly in the 
lower stratosphere, but still having substantial disagreement with observations. Including 
both orographic and non-orographic gravity wave effects, both the zonal mean zonal 
winds and zonally-averaged temperatures are closer to the ERA-40 climatology, but the 
degree of agreement between model results and ERA -40 is quite different in the 
Northern and Southern Hemispheres. We will consider this point more in the next 
section. 

4. Non-uniformity of Non-Orographic Gravity Wave Sources 

While the inclusion of orographic and non-orographic gravity wave treatments 
have brought the GISS model results closer to ERA-40 climatology results, several 
problems remain. Perhaps the most notable is seen in figure (4), where the GISS- 
OG&NOG July pole-to pole temperature gradient at 1 hPa is about 70 °C whereas in 
ERA-40, it is only about 40 °C. This is manifested in the July zonal mean zonal winds 
being too strong in GISS-OG&NOG in both hemispheres. Interestingly, the January 
GISS-OG&NOG pole-to-pole temperature gradient is actually consistent with ERA-40 
(35 °C compared to 40 °C), and this is consistent with weaker zonal mean zonal winds in 
the January GISS-OG&NOG Southern Hemisphere than are seen in ERA -40. This 
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suggests that a globally uniform non-orographic gravity wave scheme is not appropriate, 
a fact previous noted by Charron and Manzini (2002), Garcia et al. (2007), and Richter et 
al. (2010). 

Since the GISS-OG&NOG January pole-to-pole temperature gradient results 
agree well with ERA-40 in January, we have experimented with what globally uniform 
non-orographic gravity wave source function is needed to bring the July temperature 
gradient into agreement with observations. Figure (5) shows GISS-OG&NOG January 
and July zonal mean zonal wind and temperature distributions for B t = 0.004 kg m' 1 s' 2 . 
Note first that the July pole-to-pole temperature gradient at 1 hPa is about 50-55 °C, 
somewhat larger than the July ERA-40 value, but for this value of B t , the January 
temperature gradient at 1 hPa is only about 25 °C, which is about 1 5 °C less than 
observed. Consistent with these temperature distributions, the July Southern Hemisphere 
polar night jet is slightly stronger than in ERA-40, but the summer Northern Hemisphere 
easterlies are much too weak in figure (5). Interestingly though, there is evidence of 
equatorial westerly to easterly shear regions for this large value of B t that were not so 
evident in figures (1) and (3). The GISS-OG&NOG January zonal mean zonal winds in 
figure (5) are too weak in the winter lower stratosphere, but are of reasonable magnitude 
in the upper stratosphere, albeit located at too low latitudes. The January summer 
easterlies are much too weak for this large value of B t . The winter westerlies show no 
equatorward tilt, such as seen in ERA-40, in figure (5) in either January or July. 

Clearly then, some non-uniformity in non-orographic gravity wave source 
functions is required to bring GISS-OG&NOG results into agreement with observations, 
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but there should be physical justifications for this non-uniformity. Remembering that 
some of the physical sources for non-orographic gravity waves are convection, fronts, 
and spontaneous emission from jets, one can get an idea of what the nature of their 
temporal and spatial distribution might be. In a recent paper by Richter et al. (2010), they 
included explicitly computed gravity wave source functions in the Whole Atmosphere 
Community Climate Model (WACCM). Their physical sources for non-orographic 
gravity waves were convection and emission from frontal systems. Their resulting 
momentum fluxes at 100 hPa are shown in their figures (2) and (3). 

Conceptually then, we will only consider two types of physical sources for non- 
orographic gravity waves - convection and spontaneous emission from jets. This is 
different from Richter et al. (2010) in that their non-orographic gravity wave sources are 
convection and emission from fronts. Richter et al. (2010) indicate that their frontal 
source for gravity waves should often be co-located with jet sources, but there are two 
important differences. One is that the jet source should be at higher altitudes than the 
frontal sources, and the other is that jet sources should be more ubiquitous than frontal 
sources. Gong and Geller (2010) have performed a study in which they trace back 
gravity waves that they observe using high vertical-resolution radiosonde data. They find 
that in the cases they consider, they can use ray-tracing to establish that the source of 
these gravity waves are associated with jet imbalances. Furthermore, they use the linear 
model of Wang and Zhang (2010) to show that these source jet imbalance regions give 
rise to gravity waves have frequencies and wavenumbers that are consistent with the 
radiosonde observations. 
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For these reasons, we consider the non-orographic gravity wave sources to be at a 
pressure altitude of 100 hPa (Richter et al., 2010 have their non-orographic gravity wave 
source altitude at 600 hPa), and to have the spatial and temporal dependence that we 
expect from jet imbalance and convection sources. We consider both these gravity wave 
source functions to have the same spectral shape, and we “tune” the magnitudes of their 
momentum fluxes, Bt, to give zonal mean zonal winds and zonally-averaged 
temperatures that are consistent with ERA -40. The wind and temperature results are 
shown in figure (6). Figure (7) shows the temporal and spatial distribution of B t that 
were used in GISS-OG&NOG to obtain these results. This non-orographic gravity wave 
source function at 100 hPa was specified to have a Gaussian shape in latitude for both the 
tropical convective source and the extratropical jet imbalance source, and these Gaussians 
were centered at latitudes that were meant to represent the seasonal variation of the 
Intertropical Convergence Zone and the polar jet streams. The relatively broad width of 
the extratropical Gaussians are meant to represent the splitting and meandering of these 
jets. A background non-orographic gravity wave source with B t = 0.001 kg mf s' is 
taken to exist everywhere, and is superposed upon our idealizations of the tropical 
convective and jet stream sources. The background source is meant to represent any 
number of non-orographic gravity wave sources that are not associated with ITCZ 
convection and jet imbalance. This overall picture is consistent with the gravity wave 
climatologies that are derived by satellite, radiosondes and GPS data (e. g., Alexander 
and Barnet, 2007; Wang and Geller, 2003; Tsuda et al., 2000) that show greater gravity 
wave activity in the stratosphere during winter than summer. Note the time varying 
gravity wave source function shown in figure (7) is larger in the Southern Hemisphere 
than in the Northern Hemisphere . Their ratio is reasonable given the stronger Southern 
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Hemisphere storm tracks found by Guo et al. (2008) and the fact that the jet imbalance 
gravity wave source involves the square of the jet strength (see Wang and Zhang, 2010, 
for example). 

Both the GISS-OG&NOG January zonally averaged temperatures and zonal mean 
zonal winds in figure (6) agree well with ERA-40 in both the Northern and Southern 
Hemispheres in terms of the jet stream strengths and locations, the region of minimum 
shear between the tropopause jet and the polar night jet (although these July winds are 
still about 5 m s' 1 too strong), and the equatorward tilt of the polar night jet. In July, the 
wind systems agree almost as well as in January. The polar night jet shows the proper 
equatorward slope, and the region of minimum shear between the lower and upper jet 
systems in the Southern Hemisphere agrees reasonably with ERA-40. Consistent with 
the thermal wind relationship, the temperatures also show good agreement with ERA-40. 
In January, the lower stratosphere minimum is a little colder than in ERA-40 and is 
located at the same altitude. The January pole-to-pole temperature gradient is about 45 
°C in the GISS-OG&NOG model compared to 40 °C in ERA-40, and the stratopause 
temperature distributions are very similar. In July, the modeled minimum winter lower 
stratospheric temperatures compare well with ERA-40, and they occur at the same 
altitude as in ERA-40. The modeled summer stratopause temperatures are about 5 °C 
cooler than in ERA-40 in January and about 10 °C cooler in July, while the modeled 
winter stratopause temperatures are about 10 °C cooler than in ERA-40 in January and 
agree well in July. This implies that the modeled pole-to-pole temperature gradients are 
about 5 °C too high in January and agree well in July with those in ERA-40. 
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Schmidt et al. (2006) show the comparison between their M20, M23, and F20 
GISS Model-E modeled zonal mean zonal winds with CIRA (COSPAR International 
Reference Atmosphere) results in their figure (16). To examine the changes that result 
from including both orographic and non-orographic gravity waves compared to the J-drag 
used in Schmidt et al. (2006), we have preformed an identical AMIP-style run using the 
J-drag formulation of Schmidt et al. (2006). These results are shown in figure (8). 
Looking first at the zonal mean zonal winds, both the summer and winter stratospheric 
winds are much too weak in both January and July. Consistent with this, the modeled 
January and July pole-to-pole temperature gradients at 1 hPa are much too weak 
compared to observations. The winter minimum polar night temperatures are also seen to 
be too warm by about 5 °C in January and about 10 °C in July. Clearly, the J-drag is 
much too dissipative. 

It should be noted that comparing results obtained using J-drag values, that were 
“tuned” for the GISS Model-E in Schmidt et al. (2006), in the newer model for which our 
gravity wave treatment was “tuned” is not a completely fair comparison. However, it 
should also be noted that the J-drag results shown here show very similar deficiencies to 
those of the M20 and F20 results shown in Schmidt et al. (2006), so we think we are 
identifying improvements that are mainly due to using our gravity wave treatments 
instead of the J-drag to represent gravity wave effects. 

Noting that all the gravity wave “tuning” has been focused on getting reasonable 
troposphere/stratosphere zonal mean zonal wind and temperature structures, while using 
approaches that are defensible in terms of atmospheric physics, we now look at some sea- 
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level pressure results to see how they compare with the earlier results of Schmidt et al. 
(2006). In doing so, it should be noted that there have been other improvements in the 
physical treatments in the GISS model, so one cannot attribute improvements in model 
results as being strictly due to our gravity wave treatments. On the other hand, there is a 
clear indication that by getting the troposphere/stratosphere zonal mean zonal wind and 
temperature structure, we have also improved the troposphere/surface simulation. 

Figures (9) shows our model 1980-1999 January sea-level pressure (SLP) 
climatology in comparison with the 1980-1999 ERA-40 climatology, as well as a similar 
comparison using the Model-E J-drag that was described in Schmidt et al. (2006). Note 
that panel (c) of figure (9) shows maximum positive sea-level differences between our 
model and ERA-40 of less than about 7.0 hPa everywhere except for over Greenland, and 
the Arctic in the Northern Hemisphere and a small region north of Antarctica at about 
100 °E. We do see a sizable negative SLP difference (about 10 hPa) in the North 
Atlantic, and in smaller regions where there is very high topography such as in the 
Himalayas and the Andes. The J-drag differences shown in figure (9d), on the other hand, 
while being comparable at midlatitudes are much worse at high latitudes. This is likely 
due to too much polar downwelling induced by the excessive dissipation. Figure (10) 
shows a similar comparison for July. Looking at the July SLP, we see that GISS 
OG&NOG has SLP too large by up to about 12 hPa over Greenland and north of 
Antarctica at about 0 °E. Large negative differences (up to about 1 5 hPa) from ERA-40 
are confined to high topography regions over the Himalayas and Antarctica, with smaller 
departures over the Rockies and Andes regions. Quantitative comparisons of SLP in 
these regions of high topography are questionable, however, since that they depend on 
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methods for hypsometric correction. The J-drag SLP differences from ERA-40 in July 
are generally larger, again particularly at high latitudes. 

5. Toward Climate Interactive Gravity Wave Sources 

The agreement between the GISS-OG&NOG model results and ERA-40 is very 
encouraging. These results were obtained with climatological specifications of the 
gravity wave sources, but it suggests a strategy to make these specifications interactive 
with a changing climate. This is to make the jet stream spontaneous emission sources 
dependent on the model-generated jet stream strengths and locations. There are various 
alternatives for this. One could calculate nonlinear imbalance terms (e. g., Medvedev and 
Gavrilov, 1995; Plougonven and Zhang, 2007) and launch gravity waves from these 
regions. One could use the Wang and Zhang (2010) methodology to determine the 
strength of the emitted gravity waves, and use the results to construct jet-dependent 
gravity wave momentum fluxes. One can also use the Beres et al. (2005) results to make 
an interactive parameterization for the convectively generated gravity waves in the same 
way as was done in Richter et al. (2010). 

6. Concluding Comments 

In this paper, we have only compared modeled zonal mean zonal winds, zonally- 
averaged temperatures, and surface pressure distributions. We are now in the process of 
performing more extensive diagnostics on our model results. This is underway, and is 
showing promising results. This will be the subject of a companion paper. 

We want to stress that our formulation for the parameterization of unresolved 
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gravity waves has not been “tuned” for good simulations of parameters like sea-level 
pressure, or to obtain good stratospheric transports. Rather, we have implemented our 
parameterizations in a physically reasonable manner, and we have “tuned” these to get 
reasonable simulations of troposphere-stratosphere zonal mean zonal wind and 
temperature climatologies. We have taken this approach since these are the parameters 
that are directly impacted by these parameterizations. It is very encouraging that this 
approach gives good results for sea-level pressure and constituent transports (as will be 
shown in a future paper). 

There are more improvements to be implemented. For instance, we will be 
implementing a more state-of-the-art treatment of orographic gravity waves, and this may 
lead to some changes in the parameter settings for the non-orographic gravity waves. We 
will also implement climate dependent gravity wave source functions. Finally, we will 
be investigating the influence of our model top since this is known to affect the structure 
of the resolved waves and their effects on the mean flow. 
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Figure 1 . January zonal mean zonal wind from ERA-40 (a), GISS-ND (b), GISS-OG (c), GISS-OG&NOG 
(d). Wind intervals are 5 ms' 1 . Solid (dashed) lines denote westerly (easterly )winds. 
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Figure 3. Same as figure (1), but for July. 
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Figure 4. Same as figure (2), but for July. 
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Figure 5. Zonal mean zonal winds for January (a) and July (b) from GISS-OG&NOG with B t = 0.004 kg m' 1 s' 2 . The 
corresponding zonally-averaged temperatures are shown for January (c) and for July (d). 
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Figure 7. Spatial and temporal distribution of B t used for the spatially and temporally 
varying non-orographic gravity wave source function in GISS-OG&NOG to 
generate the results in figure (6). 
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Figure 8. Same as figure (6), but using the J-drag in Schmidt et al. (2006). 









Figure 9. January sea-level pressure (in hPa) distribution from GISS OG&NOG for 


1980-1999 (a); ERA-40 for 1980-1999 (b); their difference (c); and the 


difference from ERA-40 when J-drag is used (d). 
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Figure 10. Same as figure (9), but for July. 
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